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It has long since been argued that the metallic states of the single-band Hubbard Hamiltonian 
TL in two spatial dimensions (i.e. for d = 2) should be non-Fermi liquid, a possibility that would 
lead to the understanding of the observed anomalous behaviour of the doped copper-oxide-based 
superconducting compounds in their normal metallic states. Here we present a formalism which 
enables us to express, for arbitrary d, the behaviour of the momentum-distribution function n CT (fc) 
pertaining to uniform metallic ground states of TL close to S F:C r (the Fermi surface of the fermions with 
spin index a, a =|, J.) in terms of a small number of constant parameters which are bound to satisfy 
certain inequalities implied by the requirement of the stability of the ground state of the system. 
These inequalities restrict the range of variation of n a (k) for k infinitesimally inside and outside the 
Fermi sea pertaining to fermions with spin index a and consequently the range of variation in the 
zero-temperature limit of n CT (fc) for k on S F - a - On the basis of some available accurate numerical 
results for n CT (fe) pertaining to the Hubbard and the t-J Hamiltonian, we conclude that, at least in 
the strong-coupling regime, the metallic ground states of TL in d = 2 cannot be Fermi-liquid nor 
can they in general be purely Luttinger or marginal Fermi liquids. We further rigorously identify 
the pseudogap phenomenon, or 'truncation' of the Fermi surface, clearly observed in the normal 
states of underdoped copper-oxide-based superconductors, as corresponding to a line of resonance 
energies (i.e. these energies strictly do not relate to quasiparticles) located below the Fermi energy, 
with a concomitant suppression to zero of the jump in n CT (fc) over the 'truncated' parts of the 
Fermi surface. Our analyses make explicit the singular significance of the non-interacting energy 
dispersion underlying TL in determining the low-energy spectral and transport properties of the 
metallic ground states of TL. 
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§ 1. Introduction 

The unconventional behaviour of the normal states of 
the doped copper-oxide-based high-temperature super- 
conducting compounds (hereafter referred to as cuprates) 
[1,2] has necessitated reconsideration of the conventional 
theory of metals [3,4]. The latter relics on the fundamen- 
tal premise that in these the low-lying single-particle ex- 
citations are analogous to those of non-interacting metal- 
lic systems and that the electron-electron interaction 
merely gives rise to finite renormalization of the low- 
energy or low-temperature response and transport prop- 
erties of the corresponding non-interacting systems. For 
the last approximately four decades, this perspective, 
which constitutes the phenomenological theory of Lan- 
dau [5] , has a formal basis on the foundation of the the- 
ory of many-particle systems, which provides the recipe 
for the calculation of the phenomenological parameters of 
the Landau theory in terms of some material-specific pa- 
rameters (such as the concentration of the carriers) and 
the fundamental constants of nature [6,7]. 

It is generally considered that the single-band Hub- 
bard Hamiltonian TL for fermions [8-13] confined to a 
planar square lattice captures the most significant as- 
pects of the electronic correlation in the layered cuprate 
compounds [14,15,3] (for a review see [16]). In this work 
we expose some exact properties concerning the uniform 
ground states (GSs) and the corresponding low-lying 
single-particle excitations of this Hamiltonian for arbi- 
trary spatial dimensions d. We derive some novel expres- 
sions for n .(fe) corresponding to k in the close neighbour- 
hood of S F - a for any conceivable uniform metallic GS of TL 
in terms of some constrained parameters dcducible from 
£fe, /^de A (T (k;e) = ?in (T (fc) and j^de e A a (k;e) for 
k — ► iS F;cr , where A a (k;e) stands for the single-particle 
spectral function of the interacting system and \i for the 
'chemical potential' (see later). Evidently, since n (J (fc) 
itself underlies our expressions for n cr (fe) corresponding 
to k — > <S F;(T , these expressions are implicit functions of 
n CT (fc) for k — ► iS F:(J , implying a closer interdependence of 
the parameters determining the behaviour of n CT (fc) for 
k — > <S F;CT than is strictly relevant for our specific consid- 
erations in this work. The latter interdependence may 
be fruitfully incorporated in the design of the theoretical 
schemes that are aimed to the quantitative calculation of 
n CT (fc) and other correlation functions pertaining to uni- 
form metallic GSs of TL [17]. 

On the basis of the above-mentioned expressions 
for n (J (fe), from an available numerical result (ob- 
tained within the framework of the fluctuation-exchange 
(FLEX) approximation) concerning n CT (fc) corresponding 
to the Hubbard Hamiltonian in d = 2, with U/t = 8 
and n a = n s = 0.265 [18], we establish that the under- 
lying metallic state cannot be a Fermi-liquid (FL). Here 
n a stands for the number of fermions with spin index o 



per site (tr denotes the spin index complementary to a so 
that, for a —\, we have a = j), t for the nearest-neighbour 
hopping integral and U for the on-site, or intra-atomic, 
interaction energy in TL. We arrive at a similar conclusion 
by considering some available accurate numerical results 
concerning n a (k) pertaining to the t-J Hamiltonian for 
d = 2, with J/t = 0.4 and n a = n s = 0.4 [19,20], where 
J stands for the exchange integral which in the strong- 
coupling regime is formally equal to At 2 /U [21,22] (thus 
J/t « 0.4 formally corresponds to U/t w 10). 1 

Our considerations in this work are centred around a 
fictitious single-particle Green function Q a (k;e) which 
is based on the exact iV-particle GS IvE^o) (for sim- 
plicity assumed to be uniform) of the single-band Hub- 
bard Hamiltonian TL and two variational Ansatze con- 
cerning the (N ± l)-particle GS and excited states of TL, 
{|^jv±i;s)}- Here s stands for a compound index fully 
characterizing the correlated state \^m-s)- In general, 
s = (k, a) where A; is a wave-vector which in our con- 
siderations is confined to the first Brillouin zone (1BZ) 
corresponding to the lattice {Rj} underlying TL, and a 
(also a compound index) is referred to as the 'parameter 
of degeneracy' [24,25]. The necessity to invoke this addi- 
tional parameter reflects the overcompleteness of the set 
of single-particle excitations of TL (in fact of any inter- 
acting Hamiltonian) in regard to the underlying single- 
particle Hilbert space [25]. In what follows we employ 
s = to represent GSs, for example as in \^m-o), the 
M-particle GS of TL. 

We demonstrate that insofar as a limited number of 
characteristics of the interacting system are concerned, 
G<j(k; e) is equivalent with the exact single-particle Green 
function G a (k;e). The most relevant of the latter char- 
acteristics are 

(i) the GS energy -EW ; o, 

(ii) the Fermi energy e F , 

(iii) the Fermi surface S F - a and 

(iv) the GS momentum distribution function rv(fc). 

On account of these and the relative simplicity of the 
analytic structure of Ga-(k;s), we deduce a number of 
significant facts concerning the GS and the low-lying 
single-particle excitations of the uniform metallic GSs 
of TL. We deduce for instance that the Fermi surface 
<S F;CT , Ver, pertaining to one such state is a subset (not 
necessarily a proper subset) of the Fermi surface <S F °i, 



1 We should point out that the strong-coupling limit of the 
Hubbard Hamiltonian involves in addition to the t-J Hamil- 
tonian a three-site term, proportional to t 2 /U. In identify- 
ing the t-J Hamiltonian with the strong-coupling limit of the 
Hubbard Hamiltonian, one therefore discards the latter term, 
a practice which has no a priori justification away from half- 
filling. See [16] and [23]. 
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Vcr, that is that determined by the non-interacting en- 
ergy dispersion eu and the number N a of particles with 
spin index a corresponding to the interacting (N a + N s )- 
particle uniform GS of H; leaving aside the possibility 
that the set S F ^l\S F - a of k points may be non-empty (see 
later), interaction is seen to determine the topology of 
S F;a solely through determining the subdivision of N into 

N a and N s . In cases where <S F:cr is a proper subset of S F ^ , 

S F ®l\S F;r7 constitutes the pseudogap region of the puta- 
tive Fermi surface of the interacting metallic state. Our 
finding that <S F:cr C S F ?„, directly leads us to conclude 
that the Luttinger theorem [26,27] concerning the 'vol- 
ume' of the interacting Fermi sea unreservedly applies 
to the Fermi sea pertaining to uniform metallic GSs of 
H. It further leads us to uncover a kinematic constraint 
that in many instances bars partially spin-polarized uni- 
form states from being legitimate GSs of the Hubbard 
Hamiltonian. We relegate an extensive discussion of this 
kinematic constraint to a future publication; for now 
however, we only mention that the existence of this con- 
straint clarifies the conspicuous absence of partially po- 
larized ferromagnetic regions in large portions of the zero- 
temperature phase diagram of the Hubbard Hamiltonian, 
specifically in the strictly tight-binding limit where direct 
hopping of fermions is restricted to their nearest neigh- 
bours (see, for instance, [28,29]). We thus mainly ascribe 
the failure of the Stoner criterion in its incorrect predic- 
tion of uniform ferromagnetic GSs in regions of the latter 
phase diagram where there is none (see, e.g., Fig. 5 in 
[30]), to the disregard by the Stoner theory of the above- 
mentioned kinematic constraint. 

We further deduce that the FL metallic state breaks 
down at locations of the interacting Fermi surface S F[rT co- 
inciding with saddle points of the non-interacting single- 
particle energy dispersion Sk and that the two charac- 
teristic features of the pseudogap region of the putative 
Fermi surfaces of interacting metallic states are insepara- 
bly interdependent. To clarify the latter aspect, we point 
out that experimentally pseudogap regions of the puta- 
tive Fermi surfaces of the (underdoped) cuprate com- 
pounds [31-34] are characterized by: 

(a) the absence of the quasiparticle peak in the single- 
particle spectral function at the Fermi energy and 

(b) the presence of a so-called 'leading-edge' peak in 
the same function at a lower energy. 

We show that these two aspects are inviolably related. 

§ 2. Preliminaries 

For the single-band Hubbard Hamiltonian [9,12] corre- 
sponding to a rf-dimensional lattice of N L sites we have 

« = EE T «i^' + ->< Y,Y."^">- w 

i,j a i a 



where cJ. CT , Cj ;cr are canonical creation and annihilation 
operators respectively for fermions with spin index a at 
site i corresponding to site vector R i} ni :lJ :~cl. a Ci ;cr is a 
partial site number operator, U the on-site interaction 
energy (in this work assumed to be positive) and 

^£^ fe( * " ( 2 ) 

L k 

a hopping matrix element corresponding to the non- 
interacting (spin-degenerate) energy dispersion For 
simplicity, in this work we only consider the GSs of T~L 
which preserve the full translational symmetry of the un- 
derlying Bravais lattice = 1, . . . , iV L } so that both 
in Eq. (2) above and later, ^ fc signifies a sum over iV L 
points in the first Brillouin zone (1BZ) corresponding to 
the lattice {Ri}. It is common practice to choose the 
origin of the energy to coincide with T^j = N^ 1 £/-. 
By doing so, in the strict tight-binding limit, Tjj is non- 
zero only when Ri and Rj are nearest neighbours, in 
which case it is identified with —t < 0. None of our 
following considerations relies on a strict tight-binding 
approximation to {Tij}. Since the partial number oper- 
ators N a — J2i Ci [a c i: a, Vcr, commute with H, the uniform 

states of TL can be characterized by n a :=N a /N h = n—n s , 
where n:—N/N L with N = N a . 

§ 3. The single-particle Green function and some 
exact sum rules 

We seek to investigate the single-particle properties 
of the many-body Hamiltonian H as embodied by the 
single-particle Green function G a {e) [35], with e the ex- 
ternal energy parameter. For the coordinate represen- 
tation of G a (z), the analytically continued G a {e) into 
the physical Riemann sheet of the complex z plane (see 
later), we have 

G CT (r, r'; z) = £ G CT (fc; z) ^ fc (r)^(r'), (3) 
fc 

where V'fc( r ) is the normalized Bloch function corre- 
sponding to £fc and 

s z t s . a 

The physical Green function G a (k; e), with e real, is ob- 
tained through G IJ {k\e) = lim^o G a (k; e ± irf), e<^, 
where [i stands for the 'chemical potential' to be spec- 
ified below. With N a ±\,N s \s) || s} the complete set 
of (N a ± 1 + 7V ff )-particle eigenstates of Tt (the charge- 
neutral system accommodates (N a + N s ) = N particles) 
and the corresponding eigenenergies {-E/v .±i,JV ff ;s II s }, 
the Lehmann amplitudes {f s;(T (fc)} and the associated 
single-particle energies {e S :<r} are defined as follows: 
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!<*JV o .-l,JV e ;*|Cfe ;<T |*JV;0) ) £s:a < 
(5) 
(^N;0\Ck;<r\^N„ + l,N s ;s), > Mi 

{En-S) — EN a -l,N s ;s, £s;cr < M> 
(6) 
EN a +l,N g ;s - E N -fl, £s;ct > M> 

where Cfc :cr = iV L 1 zlj&jw exp(— zfc • Rj), /J, is defined 
in Eq. (8) below, and I'J/jv;!)} and E^-o are the short-hand 
notations for \^ n„,n s ;o) an d En„,n s ;0 respectively. 

Our notation in Eq. (6) is indicative of {e S :<r} consisting 
of two disjoint subsets, an aspect which has its root in 
our assumption concerning the stability of the (N a +N & )- 
particle GS of TL. By defining 

lh{;cr'- = EN-fi - E Na _i tN -. , 

^■a~EN a +l,N s -fi - E N;0 , (7) 

the disjoint nature of the above-mentioned two sets of 
excitation energies can be expressed as follows 

< M < Vcr - ( 8 ) 

For non- metallic states, fJ^. a — ^n-u ^ s nn itc, equal to the 
fundamental gap in the single-particle excitation spec- 
trum of the system. For metallic (thus microscopic) GSs, 
lipj. a — Hji-v is infinitesimally small, scaling like 1/N P 
with p > 0; however it is non-vanishing. For these states, 
the Fermi energy e F is formally equal to p^. a (in which 
case the subscript a is redundant when N a ,N 9 ^ 0; it 
is not redundant when N„ — 0) and its undue identifi- 
cation with /x (which lies in an open interval) can give 
rise to erroneous or ill-dcfincd results. To appreciate this 
aspect, consider the single-particle spectral function 

A a (k; e) := [G a {k\ e - in) - G a (k; e + in)] 

Ztti 17J.0 

= n^|f s;ff (fc)| 2 (5( e - es;ff ), (9) 

S 

for which we have 

\ /" deA a (k;e) = n a (k), (10) 

1 r°° 

-j de A a (k;e) = 1 - n CT (fc). (11) 

For some metallic states (satisfying condition (A) intro- 
duced in § 11.2.2 below), including FLs, n (T (fc) under- 
goes a discontinuous change upon transposing k from 
infinitesimally inside to infinitesimally outside the Fermi 
sea; the amount of discontinuity is by a Migdal theo- 
rem [36,27] equal to the weight of a (5-function contri- 
bution to A a (k;e) at e = e F , the Fermi energy. From 
Eq. (10) it is observed that, for this coherent contribu- 
tion to A a (k;e) to have an impact on the behaviour of 



n ff (fe) for k — > S F - a , it is required that e F = < /i. 

Similarly, from Eq. (11) we observe that there must ex- 
ist a comparable coherent contribution to A a (k;e) at an 
e inside the interval (/z, 00); for k — > iS F;CT , the required 
5-function contribution is located at e = H%. a > /x (see 
Eq. (8)). The formal mode of expressing these observa- 
tions is firstly that, the singular points of A a (k;e) coin- 
cide with the solutions of 

e k + hY> a {k-e) =e, (12) 

where X CT (fc;e) stands for the self-energy and sec- 
ondly, that this equation has one solution at (k,e) — 
(fc^ CT ,Mjv;cr) and one solution at (fe,e) = (fe+ CT , /tx^ ;cr ), 
where k~. a and k^. a are respectively infinitesimally in- 
side and outside the Fermi sea pertaining to fermions 
with spin index a in the vicinity of fc F;(T £ <S F;(T . The ne- 
cessity for the introduction of two vectors k~. a and fc+ CT 
arises from the fact that in principle (see § 3.4 in [25]), 
Eq. (12), which corresponds to a uniform GS, cannot pos- 
sess more than one solution at any given fe; in the case of 
the Hubbard Hamiltonian, where the momentum space 
is restricted to the 1BZ corresponding to the underlying 
lattice, all functions of k are periodic with the periodic- 
ity of the 1BZ, so that through Umklapp processes the 
latter principle is not of general validity (see below). 

For an arbitrary k, Eq. (12) has either no solution 
or it has real- valued solution(s) (see above) [25]. From 
the above considerations we deduce that a singular be- 
haviour (not necessarily a discontinuity) in n a (k) at 
fe = fco^ fcjj^, with fc F;(T £ <S F;(T , implies the existence of 

at least two (real-valued) solutions £fc l CT , i = 1,2, which 
do not necessarily correspond to quasiparticle excitations 
and for which ej^ < A% ;CT and e ( ^. a > /x+ ;(7 hold. The 
validity of the latter inequalities is readily verified by 
considering Eqs. (10) and (11) along the lines presented 
above. The observation of singular behaviour in n CT (fc) 
pertaining to the GS of the Hubbard Hamiltonian for 
d =- 1 in the limit of U — > 00 at k different from ±/c F;o - 
(for instance at k = ±3fc F;(J in the case corresponding to 
n a = n s = 1/4) [37,38] implies that at these k the cor- 
responding single-particle spectral function must possess 
at least two pronounced peaks, one below z%. CT and one 
above u+. a . 

Making use of the anticommutation relations for the 
canonical fermion operators c k and Cfe ;cr , it can be read- 
ily shown that 

^f* a (fc)f s;<T (fc')=S fe , fe ', fe,fc'elBZ; (13) 

S 

however 

for U^0. (14) 

fc 

The latter is a manifestation of the aforementioned over- 
completeness of the set {f s:cr (fc)||s} [25]; one, however, 
has (for s = (fe, a) see above) 
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M fc ')L (fc , Q) for tf-O. (15) 

For our later considerations it is necessary to be aware 
of the following sum rules: 

> 

^f s %(fe)f s;<T (fe') = v<{k)S kM , k,k' e 1BZ, (16) 

> 

where stands for the sum over all s for which £ S;(T <^, 
and 

The validity of the results in Eq. (16) is readily veri- 
fied by the knowledge that H commutes with the to- 
tal momentum operator P:=h J^ fe k c£ Cfe ;<T and that 

P|*JV; ) =0. 

§ 4. Fermi sea, Fermi surface, a Luttinger theorem 
and the pseudogap 

We define the Fermi sea FS<j for the fermions with spin 
index a as follows: 

FS CT :={fc || \/G a {k; n) > 0} = {k || G a {k; ft) > 0}, (18) 

which can be shown to conform with the definition given 
by Galitskii and Migdal [39] and Luttinger [27] (see 
Eqs. (6) and (94) in the latter work). Note that, in 
view of Eqs. (8) and (9), A a (k;^i) = 0, Vfe, so that 
Im[G CT (fc; fi)] = 0, Vfe. One readily verifies that the defini- 
tion in Eq. (18) identically reproduces the definition for 

the non-interacting Fermi sea FS^ :={fc || efc :cr < s F °^} 
pertaining to the GS of the non-interacting Hamiltonian 

Ho:=^efe; CT 4;ff fi fc^' ( 19 ) 

fc.cr 

where 

£k:a--=£k + e CT , (20) 
in which the constants {co-} are chosen such that 

^e(4 0) -£fe-e CT )=AT ff , ^N a = N; (21) 

k a 

in view of our later considerations, we assume that {N a } 
corresponds to the exact iV-particle uniform GS of H. 
We further define 

FS CT :=lBZ\FS a = {k \\ 1/G a (k; f£) < 0}. (22) 

With 

S F -y.={k\\l/G a {k;e F ) = 0} (23) 



the interacting Fermi surface, it follows that tS F;cr C FS^- 
(for this reason, in the following we consider k~ a and fc F;CT 
as being interchangeable); in defining S F - a it has been as- 
sumed that Mjv. ct — Mw-ct sca l es like 1/N P , where p > and 
N — > oo, taking into account the relationship e F = 
It is directly verified that the definition in Eq. (23) repro- 
duces the definition for the non-interacting Fermi surface, 
namely S F ^:={k || e fe;(T = e F °'}. We should point out that 
S F - a as defined in Eq. (23), as well as its non- interacting 

counterpart <S F °J, need not be simply-connected, as ev- 
idenced by the fact that for d = 1, for instance, nei- 
ther S F[a nor iS f "ct is simply-connected. However, even 
for cases where S F ^ and <S F:cr are not simply-connected 
sets (either jointly or separately), the number of k points 
inside FS^ and FS CT (the 'volume' of the Fermi sea) is 
well-defined. 2 This is in fact even the case for non- 
metallic GSs for which S F - a is an empty set. An aspect 
of special importance, both generally and particularly for 
our considerations in this work, is the validity of a Lut- 
tinger theorem [27,26] according to which the number of 
k points interior to FS CT is equal to The non-triviality 
of the Luttinger theorem at issue therefore rests on the 
definition for FS CT in Eq. (18) for interacting systems, for 
which no orthonormal set of single-particle wavefunctions 
can be associated with the single-particle excitations (see 
Eq. (14) above). As an essential ingredient of the proof 
of the Luttinger theorem at issue is dependent on a term- 
by-term analysis of the perturbation series for S CT (fc;e), 
it has been variously suggested that the Luttinger the- 
orem should not be of general validity. In this work we 
demonstrate that, for the uniform metallic GSs of H, 

<V C «S$. (24) 

It follows that, even in cases where S F - a C 5 F °J, the num- 
ber of fc points in the interior of FSo- is equal to that 
of FSi , establishing the validity of the Luttinger the- 
orem in the case at hand. The finding to the contrary 
by Schmalian, efc a! [40] , is a consequence of Tf. (follow- 
ing the notation used in [40], equal to — Im£fc(w = 0)), 
as calculated by these workers, not to vanish in the 
zero-temperature limit, in strict violation of the prop- 
erty A a {k]u) = ImfE^fcj/i)] = 0, Vfe, which is 
a direct consequence of the results in Eqs. (8) and (9) 
(the second expression) above; inspection of the proof of 
the Luttinger theorem at issue [27,26] makes evident that 



2 Unless the contexts imply otherwise, our statements here 
are based on the assumption that the systems under consider- 
ation are in the thermodynamic limit so that, for the purpose 
of counting the number of fc points, one can assume a continu- 
ous (and uniform) distribution of these and thus can disregard 
sets of measure zero when applicable. 
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without Im[£ CT (fc; fi)] = 0, Vfe, with \i in the open interval 
indicated in Eq. (8) above, the counting device employed 
by Luttingcr unequivocally fails. 

The result in Eq. (24) is in conformity with our ob- 
servation [17], to be briefly discussed below (see Eq. (68) 
below and the subsequent text), that n CT (fe) pertaining to 
the uniform metallic GSs of 7i is singular (not necessar- 
ily discontinuous) for all k 6 <S F ?J- In this connection, 
we note in passing that, on the basis of the arguments 
leading to the latter observation, we [17] succeed in de- 
ducing all singularities in n cr (fc) pertaining to the GS of 
the Hubbard Hamiltonian in d = 1 in the limit of infinite 
U, as calculated by Ogata and Shiba [37] and Ogata, et 
al. [38] (the latter work considers the effects of an exter- 
nally applied magnetic field) on the basis of the Bcthc 
Ansatz solution for the GS of the Hubbard Hamiltonian. 
The fact that n CT (fc) is singular for all k s S F fi, even in 
cases where S F - a C S F fl (a proper subset) implies that, 

although for k e S F ®l\S F - i(r Eq. (12) is not satisfied for 
e = e F , there must exist an energy so, with e$ < e F , 
for which this equation is satisfied (see Eqs. (10) and 
(11) above and the subsequent text). Thus the region 
S F ®l\S F;a , when non-empty, can be identified with the 
pseudogap region, which has been detected, amongst oth- 
ers [41], in the experimental photoemission spectra con- 
cerning the cuprate compounds Bi 2 Sr 2 CaCuOs+5 [31-33] 
and La2- :E Sr :E Cu04 [34] in the underdoped region (for a 
recent review of the photoemission experiments concern- 
ing cuprates see [42]). Theoretically, these observations 
have been most recently associated with the formation 
of "charge and spin gap" and "truncation of the Fermi 
surface" into disconnected 'arcs' (in d = 2) in the neigh- 
bourhoods of the saddle-points of £&, and attributed to 
"Umklapp scattering processes across the Fermi surface" 
[43]; inspired by the physics of the M-leg ladder com- 
pounds [44-50], the present authors have identified the 
corresponding metallic state as being an 'insulating spin 
liquid' (ISL) [51,52] (see also [53-58]). We relegate a 
discussion of the physical state associated with the pseu- 
dogap phenomenon to a future publication. 

Calculations of <S F;(T based on the second-order per- 
turbation expansion of T, a (k;e) in terms of the non- 
interacting Green function indicate <S F:cr % S^l 
[59-62]. By construction, in these calculations FSo- con- 
tains the same number of k points as FSi°^ (in [61], for 
instance, by choosing the second-order change in \i to be 
given by a Fermi-surface average, the number of k points 
inside the calculated FS CT is equal to that inside FSi°' ) to 
third order in U) so that within the frameworks of these 
calculations the question with regard to the validity or 
otherwise of the Luttinger theorem has not been con- 
sidered. The essential shortcoming of the perturbative 
determination of <S F;(T along the lines followed by [59-62] 
becomes evident by realizing that one of the two fun- 



damental conditions for k to be located on S F;a is that 

lm[E a (k;ek-a)} = 0, where e k . a = e k + HY, a (k; £ fe;(T ) (cf. 
Eq. (12)), a condition not taken account of in the cited 
references; this occurs because, for FLs (note this aspect) 
to second order in U, Im[E CT (A;; e)} does not feature in the 
equation for <S F;cr . Consequently, and as can also be ex- 
plicitly verified, the S F - a obtained by [59-62] does not 
satisfy one of the necessary conditions defining S F - a . The 
gravity of violating this condition can be appreciated by 
realizing the fact that an inappropriate sign (to be con- 
trasted with magnitude) assigned to Imp^fe; £fc;cr)], fol- 
lowing an incorrect subdivision of the momentum space 
into interior and exterior of the Fermi sea, amounts to the 
instability of the GS of the system under consideration. 
It follows that S F;<7 is far more tightly constrained than 
can be determined to second order in U (more generally, 
to any finite order in U). 

§ 5. Formal considerations 

We are now in a position to introduce our formalism 
that enables us to obtain exact results concerning the 
properties (i) - (iv) enumerated at the outset of this pa- 
per. To this end, we introduce the following substitu- 
tions: 

S — 1 k, f s;cr (fc) — 1 ffe ;(7 (fc), £s;a ~ ^ £fc ;CT , 
s fcGFS,, s feGFS„ 

where, in spite of suppressing the 'parameter of degener- 
acy' a associated with s (see above), FSo- and FS CT are 
the same as in the exact theory. 

Let now &^\k) and ^^\k) be two one-to-one map- 
pings according to (see Fig. 1) 

^(fcJiFS^h-lBZ, (26) 
*S >} (fc) : FSa >-» 1BZ. (27) 

These mappings, which are functions of U, are required 
to approach the following identity mappings 

*<g(fc):FS( >~FS(°), (28) 
<)(fe):FSi 0) ^FSi 0) , (29) 

for H approaching the non-interacting Hamiltonian Ho 
as defined in Eq. (19). For weakly interacting systems, 
one therefore expects ^^(fc) and ^^\k) to resem- 
ble closely the identity mappings ^^(k) and ^^(fe), 

specifically for k deep in the interiors of FSo- and FS CT 

( < ) 

respectively. We point out that <I> ^ CT (fc) do not play any 
explicit role in our following considerations. 
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1BZ 
FS 



HFS„H 



HFS„H 





FIG. 1. Schematic representation of the one-to-one mappings 
defined in Eqs. (26) and (27). 

For fj.. CT (fc) introduced in Eq. (25), we now make the 
following choice: 



V*£(*) fc G FS - fe G 1BZ, 



V*£ (*)*»<,»(»),*> fc G FSfj , fc G 1BZ. 



(30) 



It is trivially verified that f£._(fe) reproduces the result 
in Eq. (15) for U — » and in addition yields (below 
fe,fc' G 1BZ) 

E F I; CT ( fe )V( fe ') = ^.fe'' ( 31 ) 

fcelBZ 

E_ = ( 32 ) 

fcGFS^/FS„ 

in full conformity with Eqs. (13) and (16) respectively. 

( < ) 

The significance of the mappings (fc) becomes ap- 
parent by considering for a moment the possibility that 

( < ) 

(fc) are the identity mappings from the 1BZ on to it- 
self. In this case, 6 ,<-. on the right-hand side (RHS) 

*„> (fc),fc 

of Eq. (30) would be Sl k and, as a consequence, Eq. (16) 
corresponding to J^f (ST) w0lu d be reproduced by 
Eq. (32) corresponding to the summation over k G FS CT 
(fc G FS CT ) provided k,k' G FS CT (fc,fc' G FS CT ); the RHS 
of Eq. (32) would be identically vanishing otherwise. Al- 
though, for weakly-interacting systems, vf(k) and v£{k) 
are small for k G FS CT and k G FS CT respectively, this is 
not the case in general, establishing the significance of 
( < ) 

(fc) as introduced above. Note that insofar as the 
expressions in Eqs. (31) and (32) are concerned ^^(fc) 
and <&^\k) need not be any further specified than they 
be one to one over the regions indicated in Eqs. (26) and 
(27) respectively. 



§ 6. A fictitious Green function and reproduction 
OF some exact results 



Q a {k\ z):=h 



fceiBZ 



k;cr 

z — e 



(33) 



which is to be compared with G a (k; z) in Eq. (4) above. 
Making use of the expressions in Eq. (30), we have 



n CT (fe) 



l-ry(fc) 



}. (34) 



(<)— i (<) 
Here (fc) stand for the inverse of (fc); they 

exist by our assumption with regard to the one-to-one 

property of the latter mappings. Equation (34) resembles 

the equation for the Green function in the atomic limit 

[12] in which, however, n a — N a /N L is replaced by n ff (fc). 

By analogy with the exact case, we define (cf. Eq. (18) 

above) 



TS a :={k\\l/g a {k-fi) >0}, 



(35) 



and similarly for TS a . We specify // in Eq. (42) below. 

We now demonstrate that Q a (k;e) yields the charac- 
teristics (i) - (iv) enumerated above. To this end, it 
is necessary first to determine the microscopic origin of 
{fk._(fc)} and {e ( >)_i }■ In accordance with the re- 

#<< (fe);<T 

placement s — v fc in Eq. (25) above, we denote the cor- 
responding \^N <r ±i.N a -s) by \^ Na±1 ^ Ns . k }- It is readily 
verified that, by defining \^ N ±1 N _.z) according to 



I* 



I* 



1 



l*JV;0>, 



fc G FS 



fc G FSfj , 



(36) 

f£._(fe) in Eq. (30) directly follows from the defining ex- 
pression in Eq. (5) in which the 'parameter of degeneracy' 
a in s = (fc, a) is suppressed, in accordance with the sub- 
stitutions in Eq. (25). It is further readily verified that 



(37) 

in conformity with m-,s\^ m.s') — 8s,s'- I n arriving 

at the result in Eq. (37) we have made use of the fact 

/>\ _ 

that ^ct 4 (fc) are one to one so that 5 ,>% f >% = 

#„< J (fe),*„< J (fe') 

fizz'- Equations (36) concerning the (N a ± 1 + N s )- 

particle 'eigenstates' of TL give rise to the following 'eigen- 
energies': 3 



Now we define the following fictitious single-particle 
Green function G a {k] z), associated with the suppression 
of the 'parameter of degeneracy' as indicated in Eq. (25): 



3 It can be explicitly shown [17] that in general the states 
in Eq. (36) are not eigenstates of TL and thus E N ±1 N _ are 
not in general eigenenergies of TL. 
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'N a ±l,N c 



N„±l,N, 



,sl«l* 



(38) 



Let 



^±1,^:0:=. _Jpf {^iv„±i,jv ff ;fe}' ( 39 ) 

fceFSo-/ FS„ 

for which by the variational principle we have 

EN a ±l,N s -fi < ^iV„±l,iV ff ;0- (40) 

Further let (cf. Eq. (7) and sec Eq. (52) below) 



(41) 



The 'chemical potential' fi introduced in Eq. (35) has to 
satisfy (cf. Eq. (8)) 



^N-a <i 1 < V(J - 



(42) 



From Eqs. (7), (8), (40) and (41), one readily deduces 
that 



(43) 



It thus follows that, for both metallic and non- metallic 
GSs, fj, in Eq. (35) can be identified with fi. It should be 
noted that, for cases where — A% <r scales like 1/N P , 
p > (corresponding to metallic GSs, in macroscopic 
systems), Eq. (43) in principle allows for fi%. a — 
to be finite, signifying a non-metallic GS associated with 
Q a (k;z). This is, however, not a viable possibility be- 
cause, for macroscopic metallic GSs, variational energies 
for the (N a ± 1 + N s )-particle GSs deviate from their ex- 
act counterparts only by infinitesimal amounts (see also 
the considerations following Eq. (60) below). 

From the defining expression in Eq. (6), following 
Eq. (36), for the fictitious single-particle excitation ener- 
gies we have (for the poles of the fictitious Green function 
Q„{k\z) in Eq. (34), see Eqs. (46), (47) and (48) below) 



E 



N;0 



k G FS CT , 



-fe;<7 



(44) 



N .-k ~ En-,0, fc g fs ct . 



Note that, by construction Sf,. a < jjL N . a < p for fc G FSo- 

and p < < £ fe . CT f° r G FS CT (see Eqs. (39) and 

(41)). Thus, by construction (see also the considerations 
following Eq. (48) below) 



TS a — FSg-, TS a = FS a . 



(45) 



4 The subscript in £N a ±i,N s -,o is symbolic (analogous to 
in, for instance, En-,o) and does not denote k = 0. 



From Eqs. (38) and (44) we have 

Vfc G 1BZ, 



(fe);* 



'fc;er> 



(46) 



where 

r< =E 



-> 
'k;<7 



N:0 
1 



(*JV;0|4;^C fe;(r |*JV;0>, (47) 



(^N- : o\c k .,Hclj^ N . ) - E N . . (48) 



The complete disappearance of any dependence on 

( < ) 

&a (fc) of these expressions, which is directly related 
to the assumed one-to-one nature of these mappings, is 
in full conformity with the results in Eq. (45); on the one 

hand, the definitions for {k) in Eqs. (26) and (27) 
explicitly require the knowledge of the exact FS CT while, 

on the other hand, the results in Eqs. (47) and (48), 

(<) 

which have no explicit dependence on ♦fro? (fe), establish 
this knowledge to be implicitly present in the relevant 
functions e^. a and e^. a . 



§ 7. More exact results reproduced 

We are now in a position to expose the most signif- 
icant properties of Q a (k;e). To this end, let A a (k;e) 
be the single-particle spectral function corresponding to 
Ga(k;e) defined in Eq. (33), according to Eq. (9). We 
have 

AAk;e) = h{v<(k) S(e e< a ) + u>(k) 5(e - e>„)}, 

k G 1BZ. (49) 

Note in passing that in contrast with the non-interacting 
case, where A^\k;e) = hS(e — tk\a), for all k 6 1BZ, 
A a (k;e) consists of two 6 functions, approaching, how- 
ever, (fe;e) for U — > 0; since v^(k) + ^(k) = 1 (sec 
Eq. (17)) and e^. a < [i < e^. a , it follows that A a (k;s) 
in a way mimics the behaviour of the exact A a (k; e), ac- 
cording to which the particle-particle interaction gives 
rise to transfer of spectral weight from occupied to unoc- 
cupied states [63]. From the above considerations, it is 
easily verified that A a (k;e) also exactly reproduces the 
exact results in Eqs. (10) and (11), as well as 



de e A a (k; e) = Sk + Un s . 



(50) 



The fact that A a (k;e) satisfies the combination of 
Eqs. (10) and (11) (i.e. that the e integral over (— oo, oo) 
of h~ A a (k; e) is equal to unity) as well as of Eq. (50) 
implies that the first two leading-order terms in the large- 
e| asymptotic series of Q a (k;e) arc identical with those 
of G cr (fe;e) [25]. Further, for the GS total energy E N . a 
we have 
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1 If 11 
E N;0 = g ^2 £ k n ^( fc ) + 2^ / de e ^^ CT (fc;e), (51) 

from which it follows that £/v;0> the energy corresponding 
to the GS associated with Q a (k;e), is obtained through 
replacing A a (k;e) on the RHS of Eq. (51) by A a (k;e). 
By some straightforward algebra one can then deduce 
that 



£n-,o = -EW;0, 



(52) 



which shows an interesting aspect concerning the self- 
consistency of our scheme which in part relies on the 
definition for e k . a in Eq. (44), through which it depends 
on E N . . 

From the defining expressions in Eqs. (47) and (48), 
making use of the canonical anticommutation relations 
for c k . a and Cfe ;cr , we obtain 



~-k;a 



Sk + U^-, VfcelBZ, 



(53) 



where 



^(^N:o\cl. a cl,_ q ,. a ,c p/ , a/ c k ^, a \^ N , ), (54) 



PL 



(55) 



The result in Eq. (55) reflects the following exact rela- 
tionship 

+ (1 - n -( fe )) £ fe;a = £ fe + u n 9 , (56) 

which is deduced through employing Eqs. (17) and (53). 
On replacing IvPjVjo) in Eq. (54) by an uncorrelated N- 
particle state, $jv ; o), explicit calculation reveals that 



~-k\cr 



e k + Un s for fceFS^.FS 



(0) 



(57) 



where n s pertains to |$jv;o)- Note that it is due to 
the fundamental difference in the mappings defined in 
Eqs. (26), (27) and Eqs. (28), (29) that, in the 'non- 
interacting' case, e k . a is defined over FS^ (as opposed 

to the entire 1BZ), and over FS^°\ 

Having established some fundamental (and, for our 
present considerations, relevant) aspects of Q a (k;e), we 
are now capable of exposing a number of exact proper- 
ties of the metallic GSs of the Hubbard Hamiltonian. In 
line with our above considerations, we restrict ourselves 
to GSs whose spatial periodicity is that implied by the 
underlying lattice {Rj}. With 



A CT (fe):= 



n ff (fc) 
l-rv(fc) : 



r CT (fc):= 



-fc;<r 



(58) 



from Eqs. (35) and (45), making use of Eqs. (34) and 
(46), we have 



FS CT = {fc || A ff (fc) >T a (k)}. 



(59) 



Similar expressions, in terms of A CT (fc) and T a (k), are 
obtained for FS CT and <S F:cr . Following some algebra, for 
S F - a we thus obtain 

S F , a = {k || ( £p - e<J(e+ - e> a ) = 0}, (60) 

where e F = /Ujv- ct - This together with the result in 
Eq. (42), implying e5 < £?. , allow for the follow- 
ing possibility when k F - a G S F - a (recall that, accord- 

(0) v 

F;cr )■ 

= e F + A< with 
A> 



= £f 



ing to Eq. (24), <S F;cr may be a proper subset of S F[ l): 
on the one hand ef_ — e F , ef_ 

k ¥ . a \a k T . a ,u 

A< > and on the other e?. = s£, 

with A > > 0. It is readily verified that these con- 
ditions in combination with the result in Eq. (56) im- 
ply (1 - rv(fc- CT )) A< + n CT (fc+ CT ) A> = 0, which, unless 

A> = or n (7 (k F ~. rT ) = 1 and n (7 (k F . rT ) = 0, cannot be sat- 
isfied. Thus for interacting metallic GSs, the two energy 
dispersions e k . a and e k . a coincide (up to infinitesimal cor- 
rections) with £ F for k G S F:a . In other words, S F:(T in 
Eq. (60) can be defined by either of the two conditions 
e F —e k . a = and e F — e k . a = 0. From this result it follows 
that, for metallic GSs, the infimum of £ N ±1 N _. k as de- 
fined in Eq. (38) is achieved (not necessarily exclusively) 
for all k G S F:rT . 

The above observations in conjunction with Eq. (56) 
lead us to the following significant result for the uniform 
metallic GSs of H: 



e F = ek=k F . <J + U n s , Vfc F:cr G S F - a . 



(61) 



Since £ F is a constant, the anisotropy of e k together with 
the constancy of U n s with respect to fe necessitate the 
fundamental general result presented in Eq. (24). We 
point out that, on the basis of this result, for the case of 
half-filling (corresponding to n a — n s = 1/2) where by 
convention (see text following Eq. (2)) we have e k — 
for k G S F °], from Eq. (61) we deduce that 



(at half filling), 



(62) 



which constitutes a well-known rigorous theorem for d > 
1 [64,65]. The validity of the result in Eq. (62) for d = 1, 
with n a = n s = n/2 infinitcsimally different from 1/2 
(see [66]), can be made plausible by employing Eqs. (22) 
and (23) in the work by Lieb and Wu [67], which concern 
the exact solution of the Hubbard model for d — 1 (spe- 
cific to Ek = — 2t cos(fc), where k is in units of the inverse 
lattice constant) at half-filling, from which one obtains 

that (PN-a+^N-a)/ 2 = U / 2 for &l1 U and that ' for U I °) 

l^ti-cr ~ ^N-a approaches zero exponentially rapidly; for 
instance, for U/t = 1 this gap is equal to 5.0 x 10~ 3 t, 
whereas for U/t = 0.1 it is equal to 5.9 x 10 _6 i. 
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§ 8. A KINEMATIC CONSTRAINT CONCERNING PARTIALLY 
POLARIZED UNIFORM METALLIC STATES 

The independence from a of e F together with the ex- 
pression in Eq. (61) imply the following result (Fig. 2), 
applicable to uniform metallic GSs of TL corresponding 

n a - riff = — (efe p . CT - e hr . s ) when n a ,n s ^ 0; (63) 

the conditions n a ,n s ^ imply that Eq. (63) has no 
bearing on fully ferromagnetic states, such as the Na- 
gaoka state [68,69]. Leaving aside the fact that here n a 
and riff (and in consequence, fc F;(T and fc F;5 ) correspond 
to the exact {N a + Af ff ) -particle GS of TL (as opposed 
to an uncorrelated (N a + N s )-particle state), the condi- 
tion in Eq. (63) is equivalent to that within the frame- 
work of the Hartree-Fock scheme (see, e.g., Eq. (7.6) in 
[70]). This aspect is immediately understood by consid- 
ering the fact that the Hartree-Fock formalism is a varia- 
tional approach, in which the many-body Hamiltonian TL 
is exactly taken account of and the underlying approx- 
imation solely concerns the uncorrelated nature of the 
Hartree-Fock GS; in other words, Eq. (63) could have 
directly been inferred from its Hartree-Fock counterpart 
(i.e. from Eq. (7.6) in [70]). To clarify this, let us con- 
sider Eq. (63) within the framework of the Hartree-Fock 
scheme whereby its validity is undisputed. The appear- 
ance of U on the RHS of Eq. (63) has its origin in TL, 
which, as we have mentioned above, is treated exactly 
within the Hartree-Fock framework. Consequently, ac- 
counting for the neglected correlation in the latter frame- 
work, which is implicit in the exact GS I^a^o), cannot 
entail introduction of contributions to Eq. (63) with ex- 
plicit dependence on U . It follows that the difference 
between the latter expression and its exact counterpart 
must reside in n a (and consequently n s or equivalcntly, 
k F:cr and consequently k F:S ) which is an implicit function 
of U. This demonstrates that the functional form of the 
expression in Eq. (63) is universal and indeed Eq. (63) is 
valid also within the framework of the exact theory. 5 



5 It appears that the presumed first-order transition from 
paramagnetic to fully ferromagnetic uniform GSs (and vice 
versa) within the Hartree-Fock framework must have been the 
cause for the widespread disregard of the counterpart of the 
constraint in Eq. (63) in Hartree-Fock calculations; Eq. (63) 
is always satisfied for n CT = n„ = n/2, and it does not apply 
to cases where n a = n and n& = 0, or n CT = and n„ — 
n. Our careful numerical investigations [17] reveal that there 
are instances where n„ 7^ n&, satisfying n a + ria = n, solve 
Eq. (63) and moreover yield the lowest variational energy in 
comparison with the energies of both paramagnetic and fully- 
ferromagnetic uniform states. This is specifically the case for 



IBZ 




of 



FIG. 2. Illustration of the kinematic constraint in Eq. (63). 

As we have indicated above, we relegate an extensive 
discussion of the expression in Eq. (63) to a future pub- 
lication, in which we contrast the scarcity of the solu- 
tions n a ^ rig- (satisfying n a + riff = n) to Eq. (63) with 
a wide range of results that clearly indicate a similar 
scarcity of ferromagnetic GSs for the single-band Hub- 
bard Hamiltonian in particular for hyper-cubic lattices 
and strictly tight-binding Ek for various d. For complete- 
ness, for d = 1, for the tight-binding energy dispersion 
£k = — 2icos(fc) (where k is in units of the inverse lattice 
constant) Eq. (63) takes on the form 

2t 

— riff = —— [cos(7m CT ) - cos(7m 3 )] , 

which has to be solved in conjunction with the require- 
ment n a + Uff = n. Our analytic and numerical work on 
the latter transcendental equation [17] reveals that this 
equation almost entirely accounts for the fact embodied 
by a theorem due to Lieb and Mattis [71] according to 
which the total spin S of the GS of TL for d — 1 cor- 
responding to Bk — — 2icos(£;) and U 7^ 00 is equal to 
zero. Our numerical treatment of Eq. (63) for d = 2 
corresponding to 

Ek = — 2t [cos(k x ) + cos(ky)] + At' cos(k x ) cos(ky) , 

t > 0, t' > 0, at and close to the van Hove density n VH (de- 
fined as the density for which the Fermi energy of the un- 
polarized state of the non-interacting fermions coincides 
with £fc at k — (71", 0), or any of the symmetry- related k 
points) are in conformity with the findings by Honerkamp 
and Salmhofer [72,73] based on so-called temperature- 
flow renormalization-group scheme for t'/t w 0.3; for 
tf = and general n, our results conform with the Monte 
Carlo results obtained by Hirsch [74] and the findings by 
Rudin and Mattis [75]. 



instances where ey. is expressed in terms of t and t' and for 
densities n at and close to the associated van Hove densities. 
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§ 9. The zero-temperature limit of rv(fc) for k e <S F;CT 

For metallic GSs corresponding to a £ CT (fc; e) which for 
fc G iS F;(T is continuously diffcrcntiable with respect to e in 
a neighbourhood of e = e F (see condition (A) presented 
in § 11.2.2 below), we have < Z\, ¥ <T < 1, where 



Z fe :=rv(fc F ) 



(64) 



which according to a Migdal theorem [36], amounts to 
the strength of the coherent contribution to A t7 (k;s) at 
(fc,e) = (k~. a7 e F ) (compare with A a (k;^i) = 0, Vfc). In 
cases where -Zfc P . CT ^ 0, A a (k) in Eq. (58) is similar to 
n CT (fc) discontinuous for k G S Fia . In spite of these, both 
n CT (fc) and A a (k) are well-defined functions for fc F;cr G 
<S F;(T ; we have namely 

n ff (fc F;(J ):=i [rv(fc-J + n a {k+ a )] , (65) 

and similarly for A (J (fc F;(T ). Equation (65) reflects the 
fact that in the theory of Fourier integral transforma- 
tion, principal-value integrals are defined in accordance 
with the Cauchy prescription. Owing to the intrinsic 
symmetry of the Fermi function, the thermal momentum 
distribution function approaches the value in Eq. (65) for 
the temperature T approaching zero. 



§ 10. PSEUDOGAP REVISITED; DETAILED 
CHARACTERIZATION 

Before considering the behaviour of n CT (fc) for k — > 5 F;(J 
we expose two significant aspects associated with the 
case where S F - a C S F ®1 (c/. Eq. (24)). To this end, let 

fco S S F ®l\S F[a and let k$ (fcj ) denote a vector infinites- 
imally close to fco, located inside FSo- (FS CT ). With 
a continuous function of fc, on account of Eq. (61) the 
RHS of Eq. (56) is equal to e F for k — k^. In conse- 
quence of this, with Zk a ;o defined according to Eq. (64), 
from Eq. (56) we obtain 



+ (l-n,(fe+))( e >, ;(7 - £ >_ ; J. (66) 



Since fco S <S F °j\5 F;(T , by definition e< — e> must be 



> 

a finite (negative) constant. Thus, by assuming e^. a to 
be continuous functions of fc in a neighbourhood of fco, 
we observe that the RHS of Eq. (66) must be vanishing. 
It follows that 



Zk a ;a — 0, Vfco S S F ;l\S F ;a, 



(67) 

revealing a significant aspect of sj?l\S F :<j which we have 
identified as the pseudogap region of the putative Fermi 



surface of the metallic GS of H. Now, since by the 
Hcllmann-Feynman theorem 



■V (fc) 



SE 



N;0 



Seku 



(68) 



it follows that n CT (fc) must be singular (not necessarily 

discontinuous) at all fc G 5 F 2 [17]; this assertion is based 
on the observation that to any non-vanishing variation 
$£k:cr for fc G <S F °j (in the present context 5ek- a is bound 
by the requirement that 8ek- a = Seuk-.a for all elements 
of the set {u>} of the point-group operations of the un- 
derlying lattice group) corresponds a Fermi surface <S F ^' 
(whose interior contains the same number of fc points as 
S F ®1, namely N a ) that is not capable of being obtained 

from 5 f °ct through a continuous deformation; in this con- 
nection note that the functional derivative in Eq. (68) is 
constrained by the requirement that the GS of H + STi, 
with 5H:= ^<t $ e k;<r XL Cfc ; «r^fc;^ whose energy we de- 
note by En-o+SEn-o, lie in the (A r (T -|-iV 5 .)-particle Hilbcrt 
space of Tt. Thus rv(fc) is singular at fc = fc , from which 
we conclude that in the present case Eq. (12) must pos- 
sess at least two solutions at e = , i = 1,2, for 

fc = fco G <S f °ct\S f;(T , satisfying (see the paragraph follow- 
ing Eq. (24)) 



Ji) < W > ,.+ 



(2) 



(69) 



In the event Lmp^fco; e)] = for e G [^kg- a , £p], by the 
variational principle we have (see our discussions corre- 
sponding to Eq. (43)) 



'fco;c 



< £ 



(1) 

fc ;o" 



(70) 



Similarly, in the event Im[S CT (fco; e)] = for e G 
[Mtv-ctj e fe 2 ^o-] i by the variational principle we have £^ a . a > 

£ klcr- ^ or macroscopic systems it is most likely that the 
two energies in Eq. (70) (as well as the latter two en- 
ergies) are within a small difference equal. Note that, 
since by assumption fc G S^l\S F - a , Eq. (12) cannot be 
satisfied at (fc;e) = (fco;e F ). We point out that only 

when Im[S CT (fco;e)] = for e G [4^ ; a> £ fc U does the 
single-particle spectrum possess a direct gap at fc = fco, 



whose value s 



(2) 
k a \a 



— e 



(1) 
k a ;a 



does not exceed 



k a ;<J 



-fe ;o-> 



when Im[E CT (fc ; e)] = only for e in one of the inter- 
vals [£^. CT ,£ F ] and [e F ,£fc 2 Jj is the gap in question is an 
indirect gap, not exceeding e F — £^ Q . a and — e F re- 
spectively Otherwise, although there is no true gap in 
the single-particle excitation spectrum at fco G S F ^\S F - a , 
nonetheless the single-particle spectral function A a (fco; e) 
is expected to be relatively considerable in the neigh- 
bourhoods of £fc . CT and £fe . CT and suppressed in the in- 
tervening interval, rendering the designation 'pseudogap' 
in the single-particle excitation spectrum appropriately 
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descriptive. Note that, although £ k l \ a , i — 1,2, satisfy 
Eq. (12), nonetheless since in the case at hand Zk Q -,a = 
(see Eq. (67)), the spectral function A a (ko;e) does not 
possess quasiparticle peaks at e = i = 1,2, but 

resonances. Experimentally [31-34], this feature is char- 
acteristic of the under-doped cuprates for which, in the 
normal state, in certain regions of the k space that con- 
ventionally would be considered as constituting the Fermi 
surface, the expected coherent quasiparticle peak in the 
angle-resolved photo-emission spectrum is missing, and, 
further, a gap (experimentally identified as the interval 
between e F and the 'leading edge' of the broad peak in 
the angle-resolved photo-emission spectrum) in this spec- 
trum persists for T greater than the superconducting 
transition temperature T c . 

§ 11. The behaviour of n CT (fe) for k -*■ <S F;CT 

Now we proceed with considering the behaviour of 
n CT (fc) for k close to S F:rT . Equation (53) suggests that 
@k-tj should be best described as follows: 



|Cfc;<x| ~ \B(k F . a )\ \\k - fc F;CT || 7 for k -> S F 



(77) 



(71) 



This form is general and does not impose any restriction 
on the behaviour of /35 in regions where n a (fc) ^ 0. Ow- 
ing to Eq. (60) (see also the text following this equation) 
and the result in Eq. (61), for metallic states we must 
have 

6e ;C r = n s + Ck-.a, where (k-.a - for k — > S F;cr . (72) 

Employing the following series (below, fe F;cr is the point 
at which k meets S F - a in the limit k — > S F[rT ) 

£k ~ £k F;<J + a(k F:a ) ■ (k - fc F;(T ) for k -> S F:a , 
a(k F;a ):=V k s k \ k=k ^=hv F °l, (73) 

from Eqs. (53), (71) and (72) we obtain 

£ k-,a ~ A* + a(k F:a ) ■ (k - k F:a ) + U (k- a 

for k -> S F - a . (74) 

From Eqs. (53), (55), (71) and (72) we similarly obtain 

£ k:a ~ M + a(k F . a ) ■ (k - fc F;CT ) - U A a (k) C,k-.a 



for k — > S, 



F;(T7 



(75) 



where A CT (fc) is defined in Eq. (58). Thus from Eqs. (58), 
(74) and (75) we have 



rv(fc) 



-q(fc F;(T ) ■ (fc - k F . a ) - U C, k , a 
a(k F , a ) ■ (fc - fc F;CT ) - U A CT (fc) Ck-.o 

Let now (c/. Eq. (72)) 



(76) 



We consider three cases: case I, < 7 < 1; case II, 
7=1; case III, 7 > 1. The findings corresponding to 
case I equally apply to cases where, for instance, 

|Cfc;*| ~ |S(fc F;ff )ln||fc-fc F;ff ||| ||fc -fc F;CT |r for k^S F . a . 

(78) 

In fact, the results corresponding to case I also apply 
to the case corresponding to 7 = 1 when, instead of 
Eq. (77), Eq. (78) holds; in such an event, our consid- 
erations corresponding to 7 = 1 become redundant. We 
emphasize that in principle the functional form of Cfc;<r 
for fc — ► 5 F;(7 can depend on the location of fc, whether 
inside or outside FScr. Thus, for instance, the behaviour 
°f Ck,a for fc — > iS F;(7 may be governed by two powers j K 
and 7 > , with j< 7^ 7 > , for k e FS<j and fc e FS CT respec- 
tively. We shall not explicitly deal with such instances, 
however our following analyses can be readily extended 
to cover these. 

For completeness, we mention that following some al- 
gebra (involving for instance use of the equation of mo- 
tion for Ck : a in the Hcisenberg picture) we obtain 

1 r f 1 de e A„(k\ e) 

Making use of Eqs. (61), (64) and (72), from Eq. (79) we 
readily deduce 

Zkr * = ]k F f d££ [MK,*;£)-Mkt-y,£)], (80) 

which can be of experimental significance (see, e.g., [76]) 
specifically in conjunction with the well-known result, di- 
rectly deduced from Eqs. (10) and (64), namely 

Z kF , a de [A a (k F .y,e) - A CT (fc+ CT ;e)]; (81) 

although the experimentally determined A a (kf. rT ; e) may 
be too crude for the above expressions to provide reliable 
values for Z kY „ (it is, however, natural to expect that, 
in view of the smallness of \\k^. a — k~. a \\, the difference 
A a (k^. a ;e) — A a (k F ^. cr ;e) should be relatively accurate), 
nonetheless the knowledge that < Z kr . a < 1 enables one 
to enforce, through use of Eqs. (80) and (81), some degree 
of accuracy in the analyses of the measured A a (kf. a -. e). 
It is relevant to note that Eqs. (80) and (81) remain valid 

by substituting A a herein by A a (see Eq. (49)), this on 

< 

account of the continuity of e k . a in a neighbourhood of 
S F - a and of e k . a = s F for fc e <S F:cr and e k . a > jj, for 

fceiBZ. 
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11.1. Case I: < 7 < 1 



11.2. Case II: 7 = 1 



In this case, from Eq. (76) we immediately obtain 
1 



r„(fc) 



A CT (fc) 



for k — ► S F -a- . 



(82) 



Solving for the Fermi surface iS F;cr , that is solving (see 
Eq. (59) and the subsequent text) 



A„(fc) = r CT (fc), 



from the result in Eq. (82) we obtain (note that by defini- 
tion Acr(fe) > so that the negative solution A cr (fe F;cr ) = 
— 1 has to be discarded) 



Ao-(fe F;(T ) — 1, 



(84) 



which through the defining expression for A CT (fc) in 
Eq. (58), according to which 



n CT (fc) = 



Mfc) 
l + A CT (fc)' 



results in 



n<r(fcF:<j) 



(85) 



(86) 



Note the significant fact that in the case under consider- 
ation, n CT (fc) is continuous in neighbourhood of k = fe F;(T , 
so that (see Eq. (64) and our remark with regard to the 
significance of the following result) 



0. 



(87) 



We point out that for the Hubbard (at half-filling 6 or 
U/t — ► 00) (see [37,78-80]) and the Tomonaga model (see 
[81]) for d=l, both rv(±fc F;(T ) = 1/2 and Z ±k?:a = 0. 

The result in Eq. (87) implies that in the case under 
consideration, E CT (fc;e), with k G >S F;cr , is not a continu- 
ously differcntiable function of e in a neighbourhood of 
e = e F [82]. As an explicit calculation shows [82] (see 
appendix D herein), for the one-dimensional Luttingcr 
model [83,84] we in addition have that dE cr (fc; z)/dk di- 
verges for all z as k — > ^k F -„. Equations (77) and (78), 
with < 7 < 1 and < 7 < 1 respectively, are sug- 
gestive of a similar singular behaviour in VfeS CT (fc;z) as 
k — ► <S F;cr . We therefore conclude that the identifica- 
tion with "Luttingcr-liquid" by Anderson [85,86,3] of the 
metallic states of the Hubbard Hamiltonian for (specifi- 
cally) d = 2, entails that the case considered here (con- 
cerning < 7 < 1, or, in the event Eq. (78) applies, 
< 7 < 1) should encompass all uniform metallic GSs of 
the Hubbard Hamiltonian. 



6 Concerning the continuity at k = ^fc F;CT of n CT (fc) for d — 1, 
in the half- filled case, we refer the reader to [77]. 



11.2.1. General 

For the considerations in this case, it is crucial to be 
explicit with regard to the location of k (i.e. whether 
k G FScr or k e FS CT ) for k — > S F - a . This is established 
by determining the sign of the inner product h(k F - a ) ■ 
(k — k F - a ), where n(fc F;(T ) stands for the outward unit 
vector normal to S F;a at fe = k F - a , pointing from FS CT 
to FS CT . In order to simplify our subsequent analyses, 
unless we explicitly indicate otherwise, in what follows 
we assume k — fc F;(T to be collinear with n(fc F;cr ); that is, 
for sufficiently small ||fe — k F - a \\ we assume that 

k-k F , a = T\\k-k F; 4h(k F . a ), k G FS CT , FS CT . (88) 



Accordingly, we introduce 

k F . a '.=k F - a it Kfi{k F -,j*), k J, 0. 



(89) 



Following Eq. (72), in the case at hand for k — > k F:a 
we have (cf. Eq. (77)) 

Cfe;a ~ K(kf. a ) ■ (fe - k F . a ) for k G FS CT , FS CT . (90) 

The vector b a (fe) is not a smoothly- varying function of k 
as k is transposed from inside FS a through S F - a into FS CT , 
and hence our use of kf. a as the argument of the vector 
function in Eq. (90). For conciseness, in what follows we 
employ the notation 



b^Mfc* ). 



(91) 



We point out that the possibility of ||b^ll = would 
contradict the condition 7=1 considered here. 

From the asymptotic expression for r CT (fc) in Eq. (76), 
making use of Eq. (90), for k — > S F;a we obtain 



r CT (fc) 



l a -UA a (k) bf' 



fe G FS ff , FS ff , (92) 



where 



(93) 



6J:=bJ -n(fe F;<7 ). (94) 

Below we similarly employ the notation 

A^:=A CT (fc p y. (95) 

The requirement with regard to the stability of the GS of 
the system under consideration implies the simultaneous 
satisfaction of the following two conditions: 
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which guarantee that the numerator and the denomina- 
tor of the expression on the RHS of Eq. (76) are both 
positive. In arriving at the above conditions, we have 
taken into account the facts that U > 0, A CT (fc) > (see 
Eq. (58)) and a a > 0. 

Solving Eq. (83) for AJ (through equating k in A CT (fc) 
with feJ CT ), on the basis of the expression in Eq. (92) we 
obtain the following quadratic equation 



Ub* (AJ) 2 - a<J A* - (a a + Ub?) = 0. 



(97) 



In view of the conditions in Eq. (96) and the non- 
negativity of A<j(fc), of the two solutions of Eq. (97) for 



each of the vectors k = 
following are acceptable 

AJ = 1- 



fc F:(T and fc 



£1- 



Ubf 

From this and Eq. (85) we obtain 



kj. a , only the 



(98) 



ag + UbJ > 1 
a a + 2UbJ < 2' 



(99) 



It is easily seen that for a a = (signifying fc F;cr as a saddle 
point, also referred to as a van Hove point [in regard to 
the density of the non-interacting single-particle states] , 
of e fc ), n CT (fc F F (T ) = 1/2 and for U\b^\ -> oo, n a (kf. a ) -> 
1/2, implying, through Eq. (64), that Zk ¥ , a = in the 
former case and Zfc E;<T — > in the latter. 7 Later, for 
FL metallic states (encompassing those which in spite of 
sj^l\S F - a 0, share the basic characteristics of FLs) we 
deduce Ub~ = h{\v F - a —v^l) (see Eq. (Ill) below) and 
U 6+ = —h\v F[(7 , where A is a finite positive constant. 
Thus for FLs are bounded even when U — > oo, 

since for these systems || f F:cr || is finite. Evidently, this 
observation concerning U — » oo solely demonstrates the 
self-consistency of the FL theory within the framework 
of our considerations in this paper and does not rule out 
the possibility of the breakdown of a FL metallic state as 
U is increased towards large values (for instance, through 
the change of 7 = 1 into 7 = 70 < 1). 

For completeness, we mention that quantum Monte 
Carlo results obtained by Hlubina, et al. [87] , correspond- 
ing to a 16 x 16 square lattice (for R:=2t'/t = 0.94, 
U/t = 2 and n = n VH , where t' stands for the hopping in- 
tegral corresponding to next-nearest-ncighbour sites and 
n V H for the van Hove density — for definition, see above) 
"suggest that the smearing of the momentum distribu- 
tion function n(fc) does not increase close to the VH 
[van Hove] points." For the model energy dispersion 



7 We point out that, for conventional energy dispersions, 
(Xfj ^ for n CT — > 0; that is, we similarly have Zfc p <T —* in 
the low-density limit. 



£fc = k x k y and the Fermi energy e F = at the saddle- 
point k = 0, Hlubina et al. [87] reported a finite Z F within 
the framework of the i-matrix-approximation (TMA) and 
further Im[S(fc = 0; e)] ~ —ej In 2 |e| which, although un- 
conventional for d = 2, does not imply breakdown of the 
quasiparticle picture (a fact already implicit in Z F > 0). 
The renormalization-group analysis by Gonzalez, et al. 
[88,89] for a model of two-dimensional electrons (closely 
related to the Hubbard Hamiltonian, involving quadratic 
energy dispersions close to the van Hove points), estab- 
lishes a marginal-FL [90] behaviour for the single-particle 
excitations in the vicinity of the van Hove points (see also 
[91], [88,89] and [92]). The above opposing observations 
indicate that our rigorous finding relating a a = to the 
breakdown of the FL state at k = k F:a cannot be viewed 
as being a priori evident. 

Whether n cr (fc F;(T ), as presented in Eq. (65), is equal to 
1/2, depends on the equality of 1 — n a (k~ a ) with n a (k F . r7 ); 
through some algebra we obtain 



1 - n CT (fe F;(T ) 



Ubz 



2Ub a 



(100) 



which cannot be unconditionally equal to n (T (fe F h cr ), pre- 
sented in Eq. (99) above. For completeness, in contrast 
with case I considered above, in the present case and 
provided that Zk F . a ^ 0, one can similar to the case of 
FLs (to be considered below) ascribe effective mass and 
(Fermi) velocity to low-lying single-particle excitations 
of the system; however, in contrast with the case of FLs, 
these characteristics are different for the excitations in- 
side and outside the Fermi sea FS CT . Consequently, for the 
present case, one can directly generalize, for instance, the 
semi-classical theory of metallic conduction (see Chapter 
13 in [93]). By doing so, one would obtain a generalized 
Drudc formula for the conductivity tensor which along its 
principal axes would involve two distinctive masses, cor- 
responding to the single-particle excitations inside and 
outside FS CT close to 6> F:cr . In this generalized Drude for- 
mula, the quantitative difference between the masses can 
be formally fully accounted for by introducing two dis- 
tinctive relaxation times. Further, the discontinuity in 
the present case in the 'Fermi velocity' across the Fermi 
surface S F - a gives rise to a temperature-independent con- 
tribution to the conductivity, or the associated relax- 
ation times. These observations should be relevant to 
the physics of the cuprates in view of the appearance of 
two transport relaxation times in various transport coef- 
ficients of these materials (see [94-96]). 

11.2.2. Fermi-liquid metallic states 

We consider a metallic state as being a FL provided 
that the following conditions hold [82]. 

(A) S CT (fe;e), with fc G S F - a , is continuously 
diffcrentiable with respect to e in a neigh- 
bourhood of e — e F ; 
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(B) E CT (fc;e F ) is continuously differentiable 
with respect to k in a neighbourhood of fe F;cr , 

Since 5 F ?2 does not feature in these conditions, the cases 

corresponding to S^ a \S F - a ^ are also accounted for 
here. We note in passing that satisfaction of condition 
(A) is sufficient for Zy. Y a ^ 0, so that metallic states with 
Zk Y „ 7^ are not necessarily FLs. 
For FL metallic states 



Im[E CT (fc ;e )] 



Re[E ff (fe;e)] -£ CT (fc;£ F ) 



C CT (fe)| £ -£ F nin|e- e]i 



l/ 3 



as e — > e F , where C a (k) > 0, a > and |/3| > 0. On 
account of this, neglecting Im[S CT (fe; e)], the equation for 
the single-particle excitation energies, i.e. Eq. (12), yields 
real- valued solutions; these, which we denote by £fc ;CT , 
describe the mean values of the peaks in A a (k;e) for 
(fe;£) close to (fc F;(T ;£ F ). Viewed from this perspective, 
for sufficiently small ||fc — fc F;cr ||, on account of Eqs. (12), 
(53) and (61) for FLs we must have 



V fc (. 



e k + U- 



vt{k) 



A V fc (£fc + ?iE CT (fe;£fe ; , 



(101) 



where A is a positive constant which on general grounds 
we expect to be of the order of unity (see the paragraph 
following Eq. (112) below). On the other hand, for FL 
metallic states we have (in consequence of the conditions 
(A) and (B) introduced above) 



V fc £ CT (fc;£ fc;(T ) 



where vi°l is defined in Eq. (73), and 



,(0)N 



for life — fc F ., 



-0, 
(102) 



1 



J F;a 



f Vfc£fc ;cr | fc=fc 



k^(y F °l + V fe S CT (fc;£ 



*Vlfc = fc F; , 



(103) 



in which Zk ¥ cr is defined in Eq. (64) and for which we 
have the following alternative expression: 



(104) 



= [l-h dS CT (fc F;CT ;£)/d£| e=EF 



The assumption with regard to the continuous differ- 
entiability of Ek;o in a neighbourhood of <S F;(T implies that 
we have (Fig. 3) 



fe=fe? 



UK 



a a -UA+b+. (105) 



As we demonstrate below, these results are equivalent to 
(see Fig. 3) 



= V fc £<, 



a a + Ub+ =a a ~UK-K. (106) 



1BZ 




FIG. 3. Schematic forms of ej£. a and ej£ for FL metal- 
lic states. Note that (see Eqs. (105) and (106)) ej£. CT (ej^) 
smoothly goes over into ejj ( e k. ) f° r k transposed from inside 
to outside FS CT . 



To demonstrate the equivalence of the results in 
Eqs. (105) and (106), we first consider Eq. (105) from 
which we have K = — ^-t K which makes explicit that 
K and are collinear and point in opposite direc- 
tions. Making use of the expression for A+ as presented 
in Eq. (98), one readily obtains 



K 



l 

u 



(107) 



This result, which implies that 6+ = — b a — a a /U, in 
combination with Eq. (98) yields 



a; 



Ub- 



A a a a + U b a 



(108) 



From this and b a = — A+ b^ we deduce 6 J = — A a b a , 
or a a + Ub+ = a a - A" Ub~ , which is Eq. (106). 
From Eqs. (85) and (108) we obtain 



Ubz 



a a + 2Ub, 
Thus by Eq. (65) we have 



- = 1 - n CT (fe F ). 



n a{kp- : a) — 2 ' 



(109) 



(110) 



This result is significant in that it unequivocally demon- 
strates that a uniform metallic GS of the Hubbard Hamil- 
tonian whose zero-temperature limit of n CT (fc) is not equal 
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to 1/2 for k £ S F - a , cannot be a FL. Interestingly, the 
finite-temperature (T = 0.015t/fc B ) numerical calcula- 
tions based on the fluctuation-exchange, FLEX, approx- 
imation to the self-energy operator of the single-band 
Hubbard Hamiltonian on a square lattice for d — 2 with 
U/t = 8 [18] indicate n cr (fe) to deviate from 1/2 for k £ 
«S F;CT ; for n a = n s = 0.265 and k F . a = tt(0.61024 . . . , 0) 
(in units of the inverse lattice constant), n CT (fc F;(T ) ps 0.38 
(see Fig. 4 in [18]). By considering the t-J Hamiltonian 
as the strong-coupling limit of the Hubbard Hamiltonian 
(see footnote 1), in which limit the small exchange inte- 
gral J is formally equal to At 2 /U [21,22], accurate numer- 
ical results for n CT (fc) pertaining to uniform GSs of the 
t-J Hamiltonian (deduced from a twelfth-order expan- 
sion in the inverse temperature [97] and corresponding 
to T = 0.2J/fc B ) for J/t = 0.5 and 1.0 [19] and specif- 
ically J/t = 0.4 [20] {J/t ps 0.4 formally corresponds 
to U/t ps 10.0) confirm our above conclusion that the 
underlying uniform metallic GSs cannot be FL (see Ap- 
pendix for some details concerning aspects of the n CT (fc) 
associated with the t-J Hamiltonian). As our finding in 
Eq. (24) makes evident, the interpretation provided by 
the authors [19,20], namely that the Luttinger theorem 
[27,26] concerning the number of k points enclosed by the 
S F - a associated with the underlying metallic GS, breaks 
down, cannot be the appropriate interpretation. We note 
in passing that the Fermi surface as depicted in Fig. 1 of 
[20] is the locus of the k points for which rv(fc) = 1/2 
[98]. It is important to point out that, since according to 
the numerical results in [19,20] n cr (fc) ^ 1/2 for k £ 6> F:cr , 
it is excluded that = for k £ iS F:cr , so that the un- 
derlying states must fail to be FL solely on account of 
the failure of condition (B) presented above. 

In connection with the above observations, we men- 
tion that the locus of the k points for which n CT (fc) = 1/2, 
with rv(fc) the extrapolation (based on a finite mesh of k 
points) of the finite-temperature quantum-Monte-Carlo 
results by Moreo, et al. [99], pertaining to the Hubbard 
Hamiltonian on a 16 x 16 square lattice with U/t = 4 
(at T = t/(6k B ) and for n a = n g = 0.435) almost iden- 
tically coincides with S F ®1 . As is evident from our above 
considerations, this result is not sufficient (although it is 
necessary) for the underlying metallic state to be identi- 
fied as a FL. From the available numerical data we are 
not capable of deducing whether the FL metallic state 
breaks down for U/t > U c /t, with U c a finite (positive) 
constant, although in view of the above-mentioned Monte 
Carlo results [99] it is tempting to believe that U c /t>4. 
This in particular in view of the fact that, for U/t = 4 
(and T — 0.05t/k B ) the n cr (fc) as obtained within the 
framework of the FLEX approximation [100] does not 
show any noticeable asymmetry with respect to 1/2 for 
k = kf. a (see Fig. 13 in [100]), in stark contrast with the 
n CT (fc) corresponding to U/t = 8 [18], discussed above. 

From Eqs. (101), (102) and (105) and some simple al- 



gebra, we obtain 



V ^ = \(^ l + k Ub ') S \^l 



rv(fc F;(T ) 



,(0) 



(111) 



As we have indicated earlier, the three vectors in 
Eq. (Ill) arc collinear and point in the direction of the 
outward normal vector to <S F:cr at fc F:cr . Thus Eq. (Ill) 
can be identically written as v F - a = (v F °J + U b~ /h)/\, 
from which it follows that, by the requirement of the sta- 
bility of the GS of the system, A > must hold (recall 
that through Eq. (96) we have b~ > for U > 0). From 
Eq. (Ill) we deduce that 



V F :cr 



>„(0) 



A 



> 2n CT (fc- ff )-l 



(112) 



For isotropic FLs, one has v F - a = (m c /m* a ) v^l, whereby 
v F-,a < Vf-1 implies that m* > m c ; here m c stands for the 
bare electron mass and m* for the effective mass of the 
quasiparticles with spin index a. For isotropic FLs we 
thus have A = (to* /m c )n cr (fc F : cr )/(2n cr (A: F : (T ) - l). 

It is interesting to note that, for £76~ — > 0, Eq. (109) 
yields n a (k F h . a ) 1 and n (T (fc F ^ (T ) f 1. We observe 
that n a (kf. (7 ) as calculated in this Section have the 
correct limits for the diminishing strength of interac- 
tion, in conformity with the fact that non-interacting 
metallic states are indeed (ideal) FLs. It follows that 
n (T (fc F : CT )/(2n (T (fc F : (T ) - l) — > 1 for Ub~ — > so that, in 

view of the requirement that in this limit v F:a — > v F °l, 
Eq. (Ill) is seen to imply the expected result that A — > 1 
for Ub~ — ► 0; it is a priori not excluded that at least A ps 1 
for all values of Ub~ for which the metallic state under 
consideration is a FL. We note in passing that Eqs. (105) 
and (106) imply that e^. CT and e^. CT form a Dirac-type 
spectrum for k approaching <S F;(T (see Fig. 3). 



11.3. Case III: 7 > 1 



In this case, from Eqs. (74) and (75) it follows that, for 



k:cr 



£k r .„ - £k and e; 



Both of these results imply instability of the GS of the 
system, as they imply e^. a > [i for k £ FScr and e^. a < yU 
for k £ FS<j. Consequently, 7 > 1 is not permissible. 



§ 12. Concluding remarks 

In conclusion, our considerations in this work make ex- 
plicit the considerable significance of the non-interacting 
energy dispersion Sk , in terms of which the single-band 
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Hubbard Hamiltonian TL is defined, in accurately describ- 
ing the low-energy single-particle excitations of the asso- 
ciated correlated metallic systems; not only is the shape 
of iS F;(T determined by £& (leaving aside the dependence on 
, and thus on U, of e a through Eq. (21)), but also can 
the suppression to zero of the quasiparticle weight Z^.^ 
at fc F;(T 6 iS F;CT be brought about (not exclusively) by the 
choice of £& in such a way that, for a given n a , fc F:(J coin- 
cides with a saddle point of £&. Although our theoretical 
findings combined with some available numerical results 
concerning the Hubbard and the t-J Hamiltonian reveal 
that on a two-dimensional square lattice, and at least 
for U/t>8, the metallic GSs of the single-band Hubbard 
Hamiltonian are not FL, the NFL states at issue are in 
general (that is, for a rigid band energy dispersion £& de- 
scribed in terms of nearest-neighbour and next-nearest- 
neighbour hopping integrals, t and t' respectively) neither 
pure marginal-FL nor pure Luttinger-liquid (in spite of 
the assumed homogeneity of the GS, depending on £& and 
the band fillings n a and n s , the low- lying single-particle 
excitations can be strongly anisotropic in the momentum 
space); the latter states are specifically characterized by 
single-particle spectral functions which lack a coherent 
low-energy quasiparticle peak on S F;<7 and are variously 
identified with the normal metallic states of the cuprate 
superconductors. Assuming the sufficiency of the single- 
band Hubbard Hamiltonian for describing the low-lying 
single-particle excitations of the cuprate superconduc- 
tors in their normal metallic states, our findings imply 
that further considerations based on more general en- 
ergy dispersions than those expressed in terms of t and 
t' are paramount; our considerations make explicit that 
the functional form of such energy dispersions may even 
be required to depend on the band filling n, that is, a 
rigid £fc may in general not suffice. The significance of 
the t', in particular in the underdoped (in connection 
with the pseudogap phenomenon in cuprates) and over- 
doped (in connection with the occurrence of ferromag- 
nctism and strong ferromagnetic fluctuations) regions of 
the phase diagram of the single-band Hubbard Hamilto- 
nian for d = 2 has been already recognized. 
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APPENDIX: A. CONSTRAINED OPERATORS AND 
THE ASSOCIATED MOMENTUM DISTRIBUTION 
FUNCTION 

The t-J Hamiltonian is defined in terms of constrained 
(or projected) site operators {a i;(T } which are expressed 
in terms of the canonical site operators {ci ;(T } as follows 
[63,21,22] 

a i;CT :=c i;cr (l - hi- 9 ). (Al) 

With 

A L i 

and | T^) an A^-particle normalized state, with N — 
N a + N s , we define the following momentum distribu- 
tion function 

n <T (fe):=(T A r|a^ ;(T a fc;(T |T7v)- (A3) 

In this appendix we investigate the relationship between 
n CT (fc) and the 'conventional' momentum distribution 
function n a (k), defined as 

n <T (k):={'T N \^,. <r c k . a \r N ), (A4) 

by assuming |Tjv) to be a variational Ansatz for the GS 
|*jv ; o) of the t-J Hamiltonian in the subspace of the 
A^-particle Hilbert space where site double occupancy 
is excluded, that is Ci^c^ a \^ n-,o) = 0, Vz; in contrast, 
Ci;&Ci;o\Y n) 7^ 0, for some or all i. 

Making use of the commutation relation 

[ci;<7 7 ^j';(t'] _ = $i,j$cr,<7' Ci\(7i (-^-5) 

we obtain 

4>j> = c\. a tj. a + h iJ;a , (A6) 

where 
With 

n -( fe ) = ^E< T ^i4^>i T ^) e ^ (Ri " Rj) > ( A8 ) 

i,3 

and similarly for n CT (fc), from Eq. (A6) it follows that 

n CT (fc) = rv(fc) + (5n CT (fc), (A9) 

where Sn a (k) is related to {T N \hij :(7 \T N ) according to 
an expression similar to that in Eq. (A8). 

Making use of the canonical anticommutation relations 
for Ci- a and c\. a , one readily obtains the following normal- 
ordered expression: 
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(Tjv|/ii i3 - ;CT |T/v) 



- < T Af |cl ;S .cJ ;tT Ci. & C i;CT |Tat) (Jij 

+ ( T Jvl454 ( T^;ffC J ; CT |T 7V ) 



H^N\c]. s cl.„C j . > aCj. <r \? N ) 

— (Tisr |c^ct ; _ ct ^Ci^Cj-ja-c,-^ |Tjv>- (A10) 

It is directly verified that, in the event that 
Ci;gCi l( j\T n) — 0, Vi, the RHS of Eq. (A10) is, as ex- 
pected, vanishing. One further immediately observes 
that 

(Tjv | &i,i ;(7 1 Tjv) = -(TAr|c i t . (T c|;S-Ci;5'Ci; (T |Tjv) 

= -^l(TiV-2; s |^ ;< r|Tiv)| 2 < 0, (All) 

s 

the equality sign applying for any i for which 
Ci;ffCj ;0 - 1 T^r) = 0. It follows that 



5^n a (k)<0, 



(A12) 



or equivalently (see Eq. (A9)) 



The equality signs in Eqs. (A12) and (A13) apply in cases 
where c i;ff Ci ;CT |Y " N ) = 0, Vi. 

In order to gain insight into the behaviour of 5n a (k) 
and the amount that the sum on the left-hand side of 
Eq. (A13) can deviate from n CT , we now consider the case 
where |Tjv} is replaced by the normalized uncorrelated 
iV-particle GS of Wo (see Eq. (19)), which we denote by 
$jv ; o), approximating the GS |\J/ at ; o) ■ Following some 
algebraic manipulations, the details of which we shall 
present elsewhere, we obtain 



( $ W;0|4^4; CT ^C,>|$;V;0) = {$N;0 \c]. s c]. a Cj.<fC j;(T | $jV ; o) 

= ^E e " fe ' ( ^ Ri) [-^"S 0) (fe)], (A14) 



t At 



L k 



x []^2 E ni 0) (fei)ni 0) (fc 2 )nW(fc-fc 1 + fc2) 

L fci.fca 

-n|n(°)(fe)], (A15) 

where n„ (k):={^N-o\c l k . c7 Ck-a\^N-,o) is equal to the unit 

step function, unity for k e FS^ -* and zero for k <E FS^° . 
It is readily verified that for i = j both sides of Eq. (A15) 
are vanishing, exactly as in the case corresponding to the 
correlated GS I^a^o) for which we naturally also have 



Ci;g-Ci-a\^N;o) = 0, Vi (note the spin indices). Making 
use of 



1 ^ e - ife .(H i -H 3 ) =(5 . 



(A16) 



for fii 0) (fc)H$iV;o|aLafc;-l^;o} horn Eqs. (A9), (A10), 
(A14) and (A15) we obtain 

n(°)(fc) = n<°)(fc)+n A + (n 9 - 2)n ff n< r °>(fc) 

ni° ) (fc 1 )ni 0) (fc 2 )n( r °)(fc - k, + fe 2 ). (A17) 



k± ,/C2 



From the result in Eq. (A17), one readily obtains 

^^^(fc^-n.n^O, (A18) 



in conformity with the exact result in Eq. (A12) above. 
The deviation of the RHS of Eq. (A18) from zero is a 
measure for the violation of Ci- t a-Ci;a\^N;o) — 0, Vi, for 
|^at;o) replaced by |$at ; o) (see text following Eq. (A10)). 
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FIG. 4. The momentum distribution function n^(fe) (solid 
line) pertaining to a square lattice for d = 2 for k = :r(7r,7r) 
(in units of the inverse lattice constant) corresponding to 
n a = ns = 9/20, at temperature T — 0.1t/k B - The 
curve for n^(fc) + n a n^ (broken line) should be com- 
pared with curves in Fig. 1 of [19]. We have taken ac- 
count of T by replacing n^ (ft), Vfe, a, on the RHS of 
Eq. (A17) by the non-interacting Fermi function introduced 
in Eq. (A19). For completeness, here k F:a corresponds to 
%f;<j = 0.485686... (indicated by the vertical broken line); at 
Xf-.cj, fi^r (ft) + n„-n^ is approximately equal to 0.461, that is 
slightly less than 1/2. Further, in the vicinity of x = 0.436, 
nS 0) (fc) + n a n s » 0.4036 + 0.2025 = 0.6061 > 1 - n CT = 0.55, 
in violation of the inequality in Eq. (A24). 

In Fig. 4 we depict (fc) as expressed in Eq. (A17) 
above in which we have, however, replaced all n„{k), 
Vfc, <t, by the non-interacting Fermi function 
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fa(k):= 



1 



l + cxp([e k;rT -ei 0) ]/[k B T}) 



(A19) 



The function depicted in Fig. 4 concerns fc along the di- 
agonal direction of the 1BZ pertaining to a square lattice 
and corresponds to n a = n s = 9/20 and T = 0.lt/k B . 
This Figure is therefore to be compared with Fig. 1 of 
[19]. We should point out that rv(fc) in both [19] and 
[20] (corresponding to I'J/jVjo}) are subject to the normal- 
ization condition A"" 1 J2k ^a(k) = n a ( c /- Eq. (A13) and 
the subsequent remark). 

By attributing —n a n s (which in the present case is 
equal to -0.2025) on the RHS of Eq. (A18) to a uniform 
underestimation of n CT (fc) by n^(fc) over the 1BZ (see 
Eq. (A9)), n? } (fc) + n CT n ff , also depicted in Fig. 4, would 
be the function to be directly compared with the func- 
tions depicted in Fig. 1 of [19]. Indeed, the upward shift 
of n^fc) by brings the resulting function into rel- 

atively good quantitative agreement with the results in 
[19]; the results are also comparable with those as calcu- 
lated by Stephan and Horsch [101] (see Table I herein; 
for an explicit comparison of the latter results with those 
based on an high-temperature expansion, see Fig. 2 in 
[97]); for instance, for J/t = 0.4 and two holes per 20 sites 
(i.e. for n = 0.9, or n a = n s = 0.45) the latter workers 
obtained 0.54681 and 0.36247 for fc at (0, 0) and (tt, tt) (in 
units of the inverse lattice constant) respectively, to be 
compared with our upward-shifted results which for the 
indicated k points are equal to 0.583 and 0.342 respec- 
tively. We should point out that the uniform increase in 
(fc) by n a n s (in order to account approximately for 
the consequence of the prohibition of double occupancy 
in the correlated state |^jv;o)) mav at some regions of the 
1BZ give rise to violation of an exact result (see Eq. (A24) 
below); according to this result, in the present case n CT (fc) 
is bound to satisfy rv(fc) < 1 — 0.45 = 0.55. 

For completeness, n CT (fc) can be associated with a 
Green function G' a (r,r';z) similar to G a {r, r'; z) in 
Eq. (3) above whose corresponding Lehmann amplitudes 
are defined exactly as in Eq. (5), however, in terms of 
afc ;cr . It follows that for the single-particle spectral func- 
tion A' a (k; e) associated with G' a (r, r'; z) (cf. Eq. (9)) we 
have 



1 f°° 

- J de A' a {k;e) = (T N \[al. a ,a k . a } + \T N }. 



(A20) 



This is exactly the result applicable to A a (k; e), with the 
exception that, for the latter function, [aj^, Sfe;o-] + on 

the RHS of Eq. (A20) has to be replaced by [c k . a , c k . a ] + = 
1, hence we have Eqs. (10) and (11), according to which 
h^ 1 f^de A a (k;e) = 1. In the present case, where 



for uniform GSs, Eq. (A20) can be written as 



(A21) 



i r°° 

- J deA' a {k;e) = l-n 9 . 



(A22) 



This result has been earlier reported by Stephan and 
Horsch [101]. Since A' a (k;e), exactly as A a (k;e), is pos- 
itive semi definite, it follows that 

i r° 

- I de A' a {k;e) < 1 - n s for any e . (A23) 
With fj,' the counterpart of \i in Eqs. (5)-(8), we thus have 



n„(fc) 



de A' a {k-e) < 1 



(A24) 



This inequality has been invoked by Putikka, et al. [19], 
adding support to the reliability of their numerical re- 
sults. We note in passing that the maximum value taken 
by ni ' (fc) in Fig. 4 is equal to 0.4036 (corresponding to 
x = 0.436) which is indeed less than 1 - 9/20 = 0.55. As 
we have indicated above (see our considerations following 
Eq. (A19)), in the present case, where n a n s = 0.2025, a 
uniform increase by n a n„ of f^ 0) (fc) over the 1BZ gives 
rise to the violation of the inequality in Eq. (A24). □ 
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